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Abstract 

Magic numbers predicted by a 3-dimensional g-deformed harmonic oscillator with Ug(3) 
D sOq(3) symmetry are compared to experimental data for alkali metal clusters, as well as 
to theoretical predictions of jellium models, Woods-Saxon and wine bottle potentials, and 
to the classification scheme using the 3n + / pseudo quantum number. The 3-dimensional q- 
deformed harmonic oscillator correctly predicts all experimentally observed magic numbers 
up to 1500 (which is the expected limit of validity for theories based on the filling of 
electronic shells), thus indicating that Ug(3), which is a nonlinear extension of the u(3) 
symmetry of the spherical (3-dimensional isotropic) harmonic oscillator, is a good candidate 
for being the symmetry of systems of alkali metal clusters. 
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Metal clusters have been recently the subject of many investigations (see |^ for 

relevant reviews). One of the first fascinating findings in their study was the appearance of 
magic numbers 0, H ||, ^ , analogous to but different from the magic numbers ap- 



pearing in the shell structure of atomic nuclei W^. This analogy led to the early description 



of metal clusters in terms of the Nilsson-Clemenger model [|T2|, which is a simplified version 



of the Nilsson model 0, 113] of atomic nuclei, in which no spin-orbit interaction is included. 



Further theoretical investigations in terms of the jellium model [T^, |T6[ demonstrated that 
the mean field potential in the case of simple metal clusters bears great similarities to the 
Woods-Saxon potential of atomic nuclei, with a slight modification of the "wine bottle" 



type ^7j. The Woods-Saxon potential itself looks like a harmonic oscillator truncated at a 
certain energy value and flattened at the bottom. It should also be recalled that an early 
schematic explanation of the magic numbers of metallic clusters has been given in terms of 
a scheme intermediate between the level scheme of the 3-dimensional harmonic oscillator 
and the square well |I| . Again in this case the intermediate potential resembles a harmonic 
oscillator flattened at the bottom. 



On the other hand, modified versions of harmonic oscillators |]T8|, [T^ have been recently 
investigated in the novel mathematical framework of quantum algebras [3D|, which are 
nonlinear generalizations of the usual Lie algebras. The spectra of g-deformed oscillators 
increase either less rapidly (for q being a phase factor, i.e. q = e*"^ with r being real) 
or more rapidly (for q being real, i.e. q = with r being real) in comparison to the 
equidistant spectrum of the usual harmonic oscillator [^, while the corresponding (WKB- 
equivalent) potentials resemble the harmonic oscillator potential, truncated at a certain 
energy (for q being a phase factor) or not (for q being real), the deformation inflicting an 
overall widening or narrowing of the potential, depending on the value of the deformation 
parameter q. 

Very recently, a g-deformed version of the 3-dimensional harmonic oscillator has been 
constructed taking advantage of the Ug(3) D sOg(3) symmetry ||2^, ^ |28|, ^ ^ |. 
The spectrum of this 3-dimensional g-deformed harmonic oscillator has been found to 
reproduce very well the spectrum of the modified harmonic oscillator introduced by Nilsson 
P^ , without the spin-orbit interaction term. Since the Nilsson model without the spin 
orbit term is essentially the Nilsson-Clemenger model used for the description of metallic 
clusters it is worth examining if the 3-dimensional g-deformed harmonic oscillator can 
reproduce the magic numbers of simple metallic clusters. This is the subject of the present 
investigation. 



The space of the 3-dimensional g-deformed harmonic oscillator consists of the com- 
pletely symmetric irreducible representations of the quantum algebra Ug(3). In this space 
a deformed angular momentum algebra, sOg(3), can be defined PBI. The Hamiltonian of 



the 3-dimensional g-deformed harmonic oscillator is defined so that it satisfies the following 
requirements: 

a) It is an sOg(3) scalar, i.e. the energy is simultaneously measurable with the g-deformed 
angular momentum related to the algebra sOg(3) and its ^-projection. 

b) It conserves the number of bosons, in terms of which the quantum algebras Uq(3) and 
sOq(3) are realized. 

c) In the limit g — 1 it is in agreement with the Hamiltonian of the usual 3-dimensional 
harmonic oscillator. 

It has been proved that the Hamiltonian of the 3-dimensional g-deformed harmonic 
oscillator satisfying the above requirements takes the form 

^. = ^^o{[iV]g^^^-^^^^Cf|, (1) 

where is the number operator and C^^^ is the second order Casimir operator of the algebra 
sOg(3), while 

is the definition of g-numbers and g-operators. 

The energy eigenvalues of the 3-dimensional g-deformed harmonic oscillator are then 



in 



E,{n, I) = huo |Ng"+i - + 1]} , (3) 

where n is the number of vibrational quanta and / is the eigenvalue of the angular momen- 
tum, obtaining the values I = n,n — 2, . . . ,0 ot 1. 

In the limit of g — > 1 one obtains limg^iEq{n,l) = huJon, which coincides with the 
classical result. 

For small values of the deformation parameter r (where g = e'^) one can expand eq. (3) 



in powers of r obtaining |23 



Eg{n, I) = HuqU — hujQT{l{l + 1) — n{n + 1)) 
- huo^T^ + 1) - ^n{n + l)(2n + 1)) + 0{t''). (4) 
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The last expression to leading order bears great similarity to the modified harmonic 
oscillator suggested by Nilsson |TE|, 113 (with the spin-orbit term omitted) 



V = ^hup^ - fiuK{jJ- Kh^ >n), '^^'^V"X"' ^^"^ 



where 

<L^>„^^(V^). (6) 



The energy eigenvalues of Nilsson's modified harmonic oscillator are [13, 14 



Eni = hun — hufi' + 1) — -n{n + 3)^ . (7) 

It has been proved that the spectrum of the 3-dimensional g-deformed harmonic os- 
cillator closely reproduces the spectrum of the modified harmonic oscillator of Nilsson. In 
both cases the effect of the l{l + 1) term is to fiatten the bottom of the harmonic oscillator 
potential, thus making it to resemble the Woods-Saxon potential. 

The level scheme of the 3-dimensional g-deformed harmonic oscillator (for hujQ = 1 and 
r = 0.038) is given in Table 1, up to a certain energy. Each level is characterized by the 
quantum numbers n (number of vibrational quanta) and / (angular momentum). Next 
to each level its energy, the number of particles it can accommodate (which is equal to 
2(2/ + 1)) and the total number of particles up to and including this level are given. If 
the energy difference between two successive levels is larger than 0.39, it is considered as a 
gap separating two successive shells and the energy difference is reported between the two 
levels. In this way magic numbers can be easily read in the table: they are the numbers 
appearing above the gaps, written in boldface characters. 

The magic numbers provided by the 3-dimensional g-deformed harmonic oscillator in 
Table 1 are compared to available experimental data for Na clusters 0, |^, |^, ^, ^ in Table 
2 (columns 2-6). The following comments apply: 

i) Only magic numbers up to 1500 are reported, since it is known that filling of electronic 
shells is expected to occur only up to this limit Q. For large clusters beyond this point 
it is known that magic numbers can be explained by the completion of icosahedral or 
cuboctahedral shells of atoms ||^. 

ii) Up to 600 particles there is consistency among the various experiments and between 
the experimental results in one hand and our findings in the other. 

iii) Beyond 600 particles the predictions of the three experiments, which report magic 
numbers in this region, are quite different. However, the predictions of all three experiments 



are well accommodated by the present model. Magic numbers 694, 832, 1012 are supported 
by the findings of both Martin et al. 0] and Brechignac et al. 0, magic numbers 1206, 
1410 are in agreement with the experimental findings of Martin et al. [Q, magic numbers 
912, 1284 are supported by the findings of Brechignac et al., while magic numbers 676, 
1100, 1314, 1502 are in agreement with the experimental findings of Pedersen et al. 

In Table 2 the predictions of three simple theoretical models (non-deformed 3- 
dimensional harmonic oscillator (column 9), square well potential (column 8), rounded 
square well potential (intermediate between the previous two, column 7) ) are also reported 
for comparison. It is clear that the predictions of the non-deformed 3-dimensional harmonic 
oscillator are in agreement with the experimental data only up to magic number 40, while 
the other two models give correctly a few more magic numbers (58, 92, 138), although they 
already fail by predicting magic numbers at 68, 70, 106, 112, 156, which are not observed. 

It should be noticed at this point that the first few magic numbers of alkali clusters 
(up to 92) can be correctly reproduced by the assumption of the formation of shells of 
atoms instead of shells of delocalized electrons PT| , this assumption being applicable under 
conditions not favoring delocalization of the valence electrons of alkali atoms. 

Comparisons among the present results, experimental data (by Martin et al. |^ (column 
2), Pedersen et al. (column 3) and Brechignac et al. (column 4) ) and theoretical 
predictions more sophisticated than these reported in Table 2, can be made in Table 3, where 
magic numbers predicted by various jellium model calculations (columns 5-8, ^ 0, |32|), 
Woods-Saxon and wine bottle potentials (column 9, ||3^), as well as by a classification 
scheme using the 3n + / pseudo quantum number (column 10, IQ) are reported. The 
following observations can be made: 

i) All magic numbers predicted by the 3-dimensional g-deformed harmonic oscillator are 
supported by at least one experiment, with no exception. 

ii) Some of the jellium models, as well as the 3n -|- / classification scheme, predict magic 
numbers at 186, 540/542, which are not supported by experiment. Some jellium models also 
predict a magic number at 748 or 758, again without support from experiment. The Woods- 
Saxon and wine bottle potentials of Ref. predict a magic number at 68, for which no 
experimental support exists. The present scheme avoids problems at these numbers. It 
should be noticed, however, that in the cases of 186 and 542 the energy gap following them 
in the present scheme is 0.329 and 0.325 respectively (see Table 1), i.e. quite close to the 
threshold of 0.39 which we have considered as the minimum energy gap separating different 
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shells. One could therefore qualitatively remark that 186 and 542 are "built in" the present 
scheme as "secondary" (not very pronounced) magic numbers. 
The following general remarks can also be made: 

i) It is quite remarkable that the 3-dimensional g-deformed harmonic oscillator repro- 
duces the magic numbers at least as accurately as other, more sophisticated, models by 
using only one free parameter {q = e"^). Once the parameter is fixed, the whole spectrum is 
fixed and no further manipulations can be made. This can be considered as evidence that 
the 3-dimensional g-deformed harmonic oscillator owns a symmetry (the Ug(3) D sOg(3) 
symmetry) appropriate for the description of the physical systems under study. 

ii) It has been remarked 0] that if n is the number of nodes in the solution of the 
radial Schrodinger equation and / is the angular momentum quantum number, then the 
degeneracy of energy levels of the hydrogen atom characterized by the same tj, + / is due to 
the so (4) symmetry of this system, while the degeneracy of energy levels of the spherical 
harmonic oscillator (i.e. of the 3-dimensional isotropic harmonic oscillator) characterized 
by the same 2r2, + / is due to the su(3) symmetry of this system. 3n + I has been used to 
approximate the magic numbers of alkali metal clusters with some success, but no relevant 
Lie symmetry could be determined (see also P3|). In view of the present findings the 



lack of Lie symmetry related to 3?7, + / is quite clear: the symmetry of the system appears 
to be a quantum algebraic symmetry (ug(3)), which is a nonlinear extension of the Lie 
symmetry u(3). 

iii) An interesting problem is to determine a WKB-equivalent potential giving (within 
this approximation) the same spectrum as the 3-dimensional g-deformed harmonic oscilla- 



tor, using methods similar to these of Ref. [22|. The similarity between the results of the 



present model and these provided by the Woods-Saxon potential (column 9 in Table 3) 
suggests that the answer should be a harmonic oscillator potential flattened at the bottom, 
similar to the Woods-Saxon potential. If such a WKB-equivalent potential will show any 
similarity to a wine bottle shape, as several potentials used for the description of metal 
clusters do [|T5|, 0, remains to be seen. 



In summary, we have shown that the 3-dimensional g-deformed harmonic oscillator with 
Ug(3) D sOg(3) symmetry correctly predicts all experimentally observed magic numbers of 
alkali metal clusters up to 1500, which is the expected limit of validity for theories based on 
the flUing of electronic shells. This indicates that Ug(3), which is a nonlinear deformation of 
the u(3) symmetry of the spherical (3-dimensional isotropic) harmonic oscillator, is a good 
candidate for being the symmetry of systems of alkali metal clusters. 
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Table 1: Energy spectrum, Eg{n,l), of the 3-dimensional g-deformed harmonic oscillator 
(eq. (3)), for hcoo = 1 and q = with r = 0.038. Each level is characterized by n (the 
number of vibrational quanta) and / (the angular momentum). 2(2/ + 1) represents the 
number of particles each level can accommodate, while under "total" the total number 
of particles up to and including this level is given. Magic numbers, reported in boldface, 
correspond to energy gaps larger than 0.39, reported between the relevant couples of energy 
levels. 
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Table 1: 



n 


/ 


E,in,l) 


2(2Z + 1) 


total 


n 


/ 
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2(2/ + 1) 
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2 


2 


9 


5 
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22 
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11 


11 
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46 
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1 


1 


1.000 


6 
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10 


8 


12.614 


34 
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9 


3 
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14 
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2 


2 
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10 
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2 





2.243 


2 


20 


9 


1 


13.336 


6 
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3 


3 


3.023 


14 


34 


10 
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26 
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11 
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1 
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6 
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10 
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4 
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18 


58 
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54 
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4 


2 
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10 


68 


10 


2 


15.316 


10 
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4 
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2 


70 


10 
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2 
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5.116 


22 


92 


11 


7 
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30 
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0.724 
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5 


3 
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14 
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0.884 






6 


6 
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26 


132 


11 


5 
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22 


854 


5 


1 


6.238 


6 


138 


14 


14 


16.779 


58 


912 






0.860 
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6 
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18 


156 


11 


3 
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14 


926 


7 


7 
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30 


186 


12 


8 


17.410 


34 
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6 


2 


7.657 


10 


196 


13 


11 


17.490 


46 


1006 


6 
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2 


198 


11 


1 
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6 
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0.502 
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22 
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15 
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62 
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34 
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26 
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12 
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9 
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38 
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0.417 
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0.475 
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22 


1432 
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17 
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70 
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0.441 
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3 
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11 
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Table 2: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data of Martin et al. 
[4] (column 2), Bj0rnholm et al. [5] (column 3), Knight et al. [6] (column 4), Pedersen et 
al. [8] (column 5) and Brcchignac et al. [9] (column 6), concerning Na clusters. The magic 
numbers provided [11] by the (non-deformed) 3-dimensional harmonic oscillator (column 
9), the square well potential (column 8) and a rounded square well potential intermediate 
between the previous two (column 7) are also shown for comparison. See text for discussion. 
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Table 3: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data of Martin et al. 
[4] (column 2), Pedersen et al. [8] (column 3), and Brechignac et al. [9] (column 4), as well 
as to the theoretical predictions of various jellium model calculations reported by Martin 
et al. [4] (column 5), Bj0rnholm et al. [5] (column 6), Brack [2] (column 7), Bulgac and 
Lewenkopf [32] (column 8), the theoretical predictions of Woods-Saxon and wine bottle 
potentials reported by Nishioka et al. [33] (column 9), as well as to the magic numbers 
predicted by the classification scheme using the 3n + I pseudo quantum number, reported 
by Martin et al. [4] (column 10). See text for discussion. 
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